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Abstract. In this paper we introduce and analyze an algorithm for continu¬ 
ous data assimilation for a three-dimensional Brinkman-Forchheimer-extended 
Darcy (3D BFeD) model of porous media. This model is believed to be accurate 
when the flow velocity is too large for Darcy’s law to be valid, and additionally 
the porosity is not too small. The algorithm is inspired by ideas developed for 
designing finite-parameters feedback control for dissipative systems. It aims 
to obtaining improved estimates of the state of the physical system by incor¬ 
porating deterministic or noisy measurements and observations. Specifically, 
the algorithm involves a feedback control that nudges the large scales of the 
approximate solution toward those of the reference solution associated with 
the spatial measurements. In the first part of the paper, we present few re¬ 
sults of existence and uniqueness of weak and strong solutions of the 3D BFeD 
system. The second part is devoted to the setting and convergence analysis of 
the data assimilation algorithm. 
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1. Introduction 

The mathematical modeling and analysis of nonlinear flows and transport pro¬ 
cesses through a porous media is a very active field of research in mathematics and 
physics due to the challenging problems in engineering and applied sciences it cov¬ 
ers. Most models of porous media are based on Darcy’s law so 11 Darcy’s equation 
has become the model of choice for the study of the flow of fluids through porous 
solids due to the pressure gradients, so much that it has now been elevated to the 
status of a law in physics ” (see, e.g., ED- Darcy’s empirical flow model represents 
a simple linear relationship between flow rate and the pressure drop in a porous 
media 

k 

up = —- Vp, 

where Uf is the Darcy velocity, k is the permeability of the porous medium, p is 
the dynamic viscosity of the fluid, and p the pressure. Any deviation from this 
scenario is termed non-Darcy flow. Roughly speaking, this law neglects the inertia 
or the acceleration forces in the fluid when compared to the classical Navier-Stokes 
equations. Also, it assumes that in a porous medium, a large body with large surface 
area of the pores, is exposed to the fluid flow so that the viscous resistance will 
greatly exceed acceleration forces in the fluid unless turbulence sets in. However, 
there are several situations where nature deviates from Darcy’s law, for instance 
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when one deals with high velocity, molecular and ionic effects or in the presence 
of some non-Newtonian fluids phenomena. In these situations, it is imperative 
to develop a more adequate models than the classical models based on Darcy’s 
law. In practice, this can be achieved using Forchheimer equation based on an 
observation of P. Forchheimer which states that the relationship between the flow 
rate and pressure gradient is nonlinear at sufficiently high velocity and that this 
nonlinearity increases with flow rate. The Darcy-Forchheimer law states 

v P = -^ v /- 7P/I v /| 2v /: 

where 7 > 0 is the so-called Forchheimer coefficient and v/ stands for the Forch¬ 
heimer velocity and pf the density. In other words, Forchheimer law assumes that 
Darcy’s law is still valid up to an additional nonlinear term to account for the 
increased pressure drop. 


The BFeD model is then based on a Darcy-Forchheimer law. It was originally 
derived in its classical configuration (a = f3 = 0, a > 0, and b > 0, see equation 
(JT]) below) in the framework of thermal dispersion in a porous medium using the 
method of volume averaging of the velocity and temperature deviations in the 
pores (see e.g. m- A discussion of the formulation, validity and limitation of 
the BFD system can be found in [33j I21| . In this paper, we consider the following 
mathematical generalization of the BFeD model 


d t u — vAu+(u-X7)u + X7p + a |zt| 2a zt + b \u\ 2 ^u = /, 
V • u = 0,, tt | t=0 = u 0 , 


(1) 


subjected either to periodic boundary conditions, with period L , and f 2 = [ 0 , L ] 3 is 
the basic periodic domain 

( u(x + L, y, z, t ) = u(x, y + L, z,t) = u(x, y,z + L,t) = u(x, y, z, t), 

l ( 2 ) 


p{x + L, y, z, t) = p(x, y + L,z,t)= p{x, y,z + L,t)= p(x, y, z , t), 


or Dirichlet no-slip boundary conditions 


u\dQ = 0 , 


(3) 


where dfl denotes the boundary of a smooth domain 17, a > /3 > 0 are constants, 
and a and b are real numbers. The unknowns are the velocity field, u, and the 
pressure, p and / is a given forcing term. In system CD, we introduced the extra 
term b\u\ 2 ^u to model a pumping, when b < 0 , by opposition to the damping 
modeled through the term a\u\ 2a u when a > 0. Notice that in the limit case 
a = b = 0, we obtain the classical Navier-Stokes system. To our knowledge, there 
are only few mathematical results concerning this model and most of them focus 
on the case i < a < 1 and /3 = 0. The continuous dependence on Brinkman and 
Forchheimer coefficients and the convergence as v —► 0 of solutions of DBF equation 
to the solutions of 


d t u + (u ■ V) u + X7 p + a \u\ 2a u + b\u\u = /, 

are studied in [ 8 j [9] [22j [23[ [26], [29] and references therein. The long time behavior of 
solutions and the existence of global attractor to CD has been studied in [25) [32. , 34, 
l35j for very restrictive values and ranges of parameters a and /3. Also, existence, 
decay rates and some qualitative properties of weak solutions are shown in [3]. 
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Recently, we also become aware of the works pang. In 0 , the system o is shown 
to be well-posed (with b = 0 ) in the whole space for a smaller range of powers a 
than the our. Our result, with periodic boundary conditions, obviously generalizes 
their result. In [20], the system 0 is investigated. The authors show existence and 
uniqueness of solutions for all a > 1, with Dirichlet boundary conditions and regular 
enough initial data. Their argument relies on the maximal regularity estimate 
for the corresponding semi-linear stationary Stokes problem proved using some 
modification of the nonlinear localization technique. 

Next, we introduce a few commonly used function spaces. By abuse of notation, 
V will refer in the case of Dirichlet boundary conditions @, to the space {u G 
C£°(fi) : V • u = 0} and to {uis a trigonometric polynomial : V • u = 0}, in the 
case of periodic boundary conditions 0. Furthermore, we introduce 

H := closure of V in L 2 (fl) and V := closure of V in H 1 ^). 

The space H is endowed with the scalar product, (-,-) H induced by L 2 (Q). The 
Hilbert space V is equipped with the scalar product {u,v)v = u, Di u)h 

in the Dirichlet case, and {u,v)v = {u,v)h + Div)u in the periodic 

case. In particular, it is well known that V CHeH' C V', with dense inclusions 
and continuous injections, see, e.g., CM- From now on we will use the notation 
| • || p for all p > 1 instead of || • || iP(n) and || • instead of || • || £ =c (R + )L2(n)) for 
lightness of the notation. In the sequel, several inequalities will involve e,; for i = 
0,1,... originating from the application of Young’s inequality. These inequalities 
will be valid for all e, £j >0 for all i = 0, 1 ,... and we will omit mentioning this. 
Eventually, let us recall the Ladyzhenskaya, Agmon’s, and Sobolev inequalities (see, 

e.g., namnsmisi]) 

Hl 4 < «i HI 2 fluIIJj!, for all u&H 1 { fl), 

Hloo < k 2 HIJji for all u G H 2 (n), 

H | 6 < k 3 HIjji, for all ueH 1 ^), 

where ki,K 2,«3 > 0 , denote dimensionless scale invariant constants that depend 
only on the shape of the domain Q. These inequalities will be used tacitly in the 
estimates we will establish in this paper. 

Our starting point is the existence and uniqueness of weak solutions to system 
0 , 0 . More precisely, we have the following 

Theorem 1.1. Let f G L°°(R + ;H) and uq G H. Assume a > (3 > 0, a > 0 and 
6g1, then systems 0.0 and 0,0 have a weak solution satisfying 

u G C 0 (K + ;H weak ) O L™ C (R+; H) nLf oc (R+;V) O Lf 0 a + 2 (R + ; T 2 q+ 2 (D)), 
and 

lim sup l'u(t) || 2 < po, 

t —^ + 00 

where 

1 
2 


Po = 



2(1 + l{b< 0 }) 


Hko} 


'4' 


m ir 

^1 , A 


+ 2 


> n 

a 


.a. 

J . 
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and l{b<o} stands for the characteristic function of the set {b < 0}. Moreover, 
if 2a > 3, then the weak solutions depend continuously on the initial data, in 
particular they are unique. 


This Theorem ensures the existence and uniqueness of weak solutions to system CD 
with Dirichlet type boundary condition d3j) or periodic boundary conditions 0 if 
a > | . This result was already obtained in nng. It is rather easy to extend 
the existence to the range a > — i but not the uniqueness. Furthermore, one can 
also show the so-called finite-time extinction of the solutions if —5 < a < |, using 
Gagliardo-Nirenberg-Sobolev inequality and we refer the reader to [3]. Let us also 
mention that the assumption / £ L°°(K. + ;H) can be weakened by assuming only 
the local integrability of time of ||/|| H . 

For more regular solutions of system 0^0 we have: 


Theorem 1.2. Let f £ L°°(R+;H) and uq £ V. Assume a > 1,0 < /3 < a, a > 0 
and b £ M. Then system CD-0 has a unique global strong solution satisfying 


u £ Cg(R+; V) nLf oc (K+;ff 2 (H) n V) nL^ c +2 (R+;L 2Q+2 (fi)), 


and limsup^.^ ||Vu(t )|| 2 < pi, where 

'1 
v 


Pi '■= 


with 


and 


Vo = 


Ai—2 


2 + (Ai + 1 ) 

u a (l + 2a) 


1 

1 + - 
v 


2(1 + l{h<o}) 


46 < 0 } 


+ {Ai + 1 ) 


4^[|&|(1 + 2/3)] c 


r , 

Vo H- 

11 

L v 

JJ 


2(1 + 1 { 5 < 0 }) 


+ l{b<o> \b\ 


m 


^ - 8 + i 


[a(l + 2a)]^ 


|fl|, rp := p 0 + 
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Furthermore, the solution depends continuously on the initial data. If in addition, 
ifu 0 £ L 2 a+ 2 {n) then u £ L 1 “(R+; L 2a+2 (tt)) and d t u £ L 2 oc {R+■ L 2 (Cl)). 


In particular, Theorem 11.21 ensures the global well-posedness of the traditional 3D 
BFeD system with periodic boundary conditions. Let us mention that we use here 
periodic boundary conditions instead of the Dirichlet ones to avoid technicalities 
and present a clear picture of the data assimilation algorithm (see below). However, 
our results remain valid with Dirichlet boundary condition as it will be done in a 
follow-up paper. In addition, we observe that these strong solutions are unique 
within the class of weak solutions as it will be described below. Observe, however, 
that as a consequence of Theorem ll.ll it follows that system JT]),0| has an absorbing 
ball in L 2 (H), and by the regularity properties guaranteed by Theorem 11.21 (for 
a > 1 ) one is able to show that our system has a finite-dimensional global attractor, 
provided f £ H is time independent (cf. [TOl l3l] ). 


Remark 1.3. Let us mention that in the paper m, there is no result concerning 
global existence and uniqueness of strong solutions with initial data in H 1 . However, 
using the maximal regularity technique, the authors showed global existence and 
uniqueness for initial data in u 0 £ H 2 . This implies in turn that the solution 
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belongs to L°°(R + ; L°°{Q)) which yields simpler estimates concerning the nonlinear 
terms. 

Theorem II.21 and Theorem 11.11 show that the nonlinear damping term, when a > 0, 
dominates the pumping term for all & 6 R. However, when a < 0, we can show that 
the pumping term may dominate in some particular situations, for large enough 
initial data, a finite-time blow-up of the solutions occurs. In other words, we show 
that it is possible to construct a finite-time blowing up solutions when a < 0 
regardless the value of b. For that purpose, we consider system © in the periodic 
channel 

& = {{x,y,z) ■ ifT, y £ T, 0 <z<L}, (4) 

for a given L > 0, subject to the mixed periodic and Diriclrlet boundary conditions 

{ u(x , y , z, t) = u{x + L, y , z, t) = u(x, y + L,z,t) = u(x + L,y + L, z,t), 

(5) 

u(x,y,0,t)=u(x,y,L,t) = 0, for all x,y £ T and z€E[0,L]. 

Then, we have the following: 

Theorem 1.4. Let be as in ©, a < 0 and b £ R. Assume that a > j3 > 0. 
Then, there exists an initial data uq for which the corresponding solution of system 
©■© blows up in finite time. 

Basically, the idea of data assimilation is to incorporate spatially discrete measure¬ 
ments and observations into the physical or numerical model. The main motivation 
of data assimilation is that deterministic models are imperfect, most of the time 
incomplete because of the inability to account for all relevant process. In particu¬ 
lar, it has played a central role in the improvement of weather forecast and there 
is a growing interest in applying data assimilation to problems arising from sev¬ 
eral applications. At the mathematical level and in the case of our model ©, the 
method of continuous data assimilation was introduced first in [4] in the context of 
the 2D Navier-Stokes equations. A follow-up work was done in ,2] and [T2] (see also 
references therein, in particular [ 6 ] for the case of stochastically noisy data). The 
algorithm can be described as follows. Assume the initial data uq of system (JT]) , d21) 
is missing. Also, assume that discret spatial observations of u(t) can be used to 
construct an interpolation operator I}, ( u ) where h denotes spatial resolution of the 
collected measurements. For instance, Ih{u) can represent an interpolant operator 
based on spatial observations of the reference solution of © at a coarse spatial 
resolution of size h. Now, Let /i > 0 be a nudging parameter, to be fixed later on, 
and consider the following data assimilation algorithm 

( d t v — v Av + (v ■ V)v + V q + a\v\ 2a v + b\v\ 20 v = / + p, (lh(u) - T h [v )), 

[ V • v = 0, u| t=0 = v 0 , 

( 6 ) 

where vq is arbitrary (to be specified in the statement of the Theorems) with the 
periodic boundary conditions 

{ v(x + L, y, z, t) = v(x, y + L,z,t) = v(x, y,z + L,t) = v(x, y, z, t), 

( 7 ) 

q(x + L, y, z, t) = q(x, y + L,z,t)= p{x, y,z + L,t)= q{x, y, z, t ), 
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Then, we can obtain the original unknown solution u(t) of system (JT]), asymp¬ 
totically in time, as the limit of v(t), the solution of ©, provided the nudging 
parameter fj, is large enough and the resolution parameter h of the measurements 
of u(t) is small enough. Furthermore, the convergence holds at an exponential rate. 

In this paper, we will consider the following types of interpolant operators. The 
first is a linear operator that approximates identity, X\ : H 1 (fl) —> L 2 {Ll), such 
that 

( 8 ) 

for every ip £ i7 1 (il), where Cq > 0 is a dimensionless constant. The second inter¬ 
polant operator, I 2 : H 2 (Q) —> L 2 (fl), is also a linear operator that approximates 
identity such that 

U - llml < Co h 2 IIV iPWl + Cl h 4 IIA m (9) 

for every ip £ H 2 (Q ), where Co > 0 and C\ > 0 are dimensionless constants. 

Before we state our results concerning system let us give examples of such 

interpolant operators. For instance, let <p k {x) = jp x , for k £ 1?. The 
projector onto the low Fourier modes, up to wave numbers |fc| such that \k\ < L 2 ^ 77 J 
defined as 

2ftWOW = p kip = ^(x), 

l fc l< L-2W7T J 

where ip(x ) = ^ fcgz3 ip k (pk( x )i is an interpolant of the first type. We leave the 
exercise of showing that this interpolant satisfies the functional inequality \ip — 
I ^(?/>)| 2 < coh 2 \\Vip\\^ to the reader. The volume element operator is also an 
interpolant of the first kind. More precisely, if we divide in cubes f2*,, for all 
k = 1,..., N, with edge LN ~3 so that |fifc| = iV -1 L 3 , then 

T hW(x) = ^ y) d Vi h = LN ~i 

k=1 WC Jn k 

That is, we consider that the average of the function ip on each cube Q k is given, 
and we refer to m for a proof of (© when ip £ H 1 (n). Eventually, the interpolant 
operator obtained by using measurement at a discrete set of nodal points in f 2 = 
[0,L] 3 . Let Xi £ f li, where Lli denotes the i th cube introduce above, be the points 
where the measurement of the velocity of the flow is taken. Therefore, the nodal 
point interpolant operator can be defined as 

N 

o) = 1 {w}( a; )- 

k= 1 

We refer the reader to ©la for a proof showing that I 2 satisfies ||ip — (VO II 2 < 

c 0 h 2 \\V ip\\l + ci h 4 | A ip\\l. 

Now, we turn to the statement of our results about system Concerning the 

first interpolant, we have 

Theorem 1.5. Let f £ L°°(R + ; H) and vo £ V. Assume a > 1, 0 < /3 < a, a > 0, 
b £ R, and 2p,co h 2 < v. Then system 0 © with interpolant has a unique global 
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strong solution satisfying 

V G Cg(R + ; V) n l£ c (R+; H 2 {M) n V) n Lf“ +2 (R + ; L 2 a+ 2 (n)). 


Moreover, ifv o G L 2a+2 (Ll), thenv G Lgc(R + ; L 2a+2 (Ll)) andd t v G L 2 oc (R + ; L 2 (fl)). 
Furthermore, the solution depends continuously on the initial data Vo in V norm. 


This Theorem ensures the well-posdness of the data assimilation algorithm when 
an interpolant of the first kind is used for all relaxation (nudging) parameter p > 0 
provided that h is small enough. Furthermore, if p is large enough, then we have 


Theorem 1.6. Let f G L°°(R + ;H), Vo G V, a > 1,0 < /? < a, a > 0, and b G R. 
Let p large enough such that 


p> 2 


2V 


+ K 4 + 1 


{ 6 < 0 } 


( 2 «o \b\) c 
{an o)P 


where K is given in Proposition \Z.1\ below. Let h be small enough such that 
2pcoh 2 < v, and u(t) be the global unique strong solution of system JT]),© given 
by Theorem U.lA Then the global unique strong solution of system with in¬ 

terpolant Ih = Tgiven by Theorem 1 1. ,51 sa.tish.es ||it(t) — v(t )|| 2 —> 0, as t —»■ +oo, at 
an exponential rate. Furthermore, ifuo,vo G V satisfying |uo|fl-i < K, ||iio|| ff i < K 
with K > y/+ p 2 , 1 < a < 2, and p is large enough so that Gl or El holds if 
1 < a < | or | < a < 2, respectively, then ||u(t) — V)lliLi(n) ~t 0, as t —> +oo, at 
an exponential rate 


Theorem E© shows that our algorithm converges toward the solution of the original 
system ©,© in L 2 norm, Moreover, Theorem E© shows the convergence in the 
Ft 1 norm if 1 < a < 2 for restricted initial data vq to be inside the absorbing ball 
of the solution of the original system. Notice that the restriction on the initial data 
Vo does not limit the applicability of our algorithm. Indeed, K depends on the size 
of the absorbing ball of u(t), \Jp\ + p 2 , which is given a priori. This restriction is 
a consequence of the fact that we are enable to control uniformly with respect to p 
the H 1 norm of v(t) for all initial data vo- When, the second interpolant is used, 
we have 


Theorem 1.7. Let f G i°°(R + ;H) and let uo,vq G V such that ||fo||jyi < K, 
ll u o||tfi < K, with K > \JPo + p 2 ■ Assume 1 < a <2,0 < /? < a, a > 0, and b G R. 
Let u(t) be the global unique strong solution of system ©,© given by Theorem \1.2\ 
Let p be large enough such that (l83l) holds if 1 < a < |, or m holds if | < a < 2, 
and h be small enough such that 15/rmax{co, yfcf} h 2 < v. Then system with 

interpolant operator Ih = I ^ has a global unique strong solution satisfying 

v G C°(R+; V) n L 2 oc (R+-H 2 (n) n V) n L^ c (R+ ; L 2a+2 m, 
and 

d t vGL 2 oc (R+-,L 2 (n)). 

Furthermore, the solution depends continuously on the initial data and satisfies 
IK*) - f(*)|| i? i (Q) —> 0, as t —> +oo, at an exponential rate. 

Observe that in this Theorem, we do not assume the initial data no to be in V 0 
L 2a+2 {VL) as in Theorem ll.51 This is because we restrict the range of a to 1 < a < 2, 
so that by Sobolev embedding vo G V implies that L 2a+2 (fl). 
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2. Proof of Theorems 11.11 and 11.21 

In this section, we prove Theorem 11.11 and Theorem 11.21 The global (in time) 
existence of solutions is shown in a classical way by proceeding in three steps. First, 
we use a Faedo-Galerkin approximation procedure to show the short time existence 
of solutions. The reader is referred to, for example, ma cEiGni for details. Next, 
we obtain the necessarya priori bounds that allow to extend the solution of the 
Faedo-Galerkin system globally in time. Eventually, we pass to the limit in the 
approximation procedure using the Aubin compactness Theorem, relying on the 
established a priori bounds. 

2.1. Existence of local solutions: Let Wk(x) be the eigenfunctions of the Stokes 
operator, corresponding to the eigenvalues A k, for k = 1 , 2 ,..., repeated according 
to their multiplicities such that 0 < Ai < A 2 < • • • < A„ < An+i < • • • and 
V m := Span{Wi,... ,W m }. It is well known that the family {Wk\kLi forms an 
orthonormal basis of H (see [10], IT (30)(31]). Now, let T > 0, be an arbitrary time, 
and consider the Faedo-Galerkin approximate function 


u m(x,t) = '^g l rn (t)W i {x) £\ m . (10) 

k =1 


Using system ©,©, the unknown coefficients g l m = J Q u m -Wi dx , for i = 1,2 ,..., m 
solve the following system of ordinary differential equations in V m 


d 

dt 


/ u„, • Wk dx + v / Vu m • VWfc dx+ (u m • V)u m • W V- dx 
J £1 J £1 Q 

+ a f |u m | 2a u m • Wfc dx + b f |u m | 2/3 u m • W k dx = [ f-W k dx, (11) 
J n Jq Jq 

ffm(t = 0) = j uo-Wfcdx, (12) 


for all k = 1,..., m. Since the vector field in system DEMED is locally Lipscliitz, 
the system admits a unique solution g!^(t) £ C' 1 ([0, T m ]), for some small interval 
[0, T m \ C [0, T], 

The solutions of system mm are obtained as the limit of a subsequence u mj (as 
j —> 00 ) of the sequence of solutions u m of the Faedo-Galerkin system mm- 
Specifically, this will be achieved by the means of the Aubin compactness Theorem, 
and the corresponding bounds are reached by Banach-Alaoglu Theorem. Also, let 
us mention that for clarity of the presentation, we will not introduce explicitly 
the Leray projector and the Stokes operator in dm, which is usually used in the 
analysis of the Navier-Stokes equation to get rid of the pressure. We will also omit 
the exercise of passing to the limit to the reader in (11011121) using the a priori bounds. 

2.2. Energy Estimates: In this paragraph, we show that u{t) is uniformly bounded 
in L 2 (Q.) with respect to time, Vit £ L 2 (0, T; L' 2 (il)) and u £ L 2a+2 (0, T; L 2 (Vl)). 
Let us mention that since the average of u is not zero, we cannot rely on Poincare 
inequality, as it is usually used in the theory of Navier-Stokes. Instead, we rely 
strongly on the damping nonlinear term as it will be clear in the sequel. For this 
purpose, we multiply the first equation of © by u and integrate over Q. Integrating 
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by parts, we obtain 


( 13 ) 


Since 0 < /? < a, we use Holder and Young inequalities to get 


Ml 2 < e MI 


2a+2 

2a+2 


HZXl<e\\u\llXl 


f ■ udx < e ||u|| 


2a+2 

2a+2 


|H|, 




2“ + l 14 

ll/llitx 


M- 


<^\\u\\Zll+\\f\\^ 

Plugging (US) and (USD in m and optimizing in e leads to 

(ll/ll2 + J 


j t 1 ^ 1 2 + 2v 1 v u \l + a mTcXI < 2 


where we set 


Vo ■= 


2(1 + l{b<o}) 


+ 1 { 6 < 0 } H 


4|6| 


|H|. 


Thanks to 0, we get (for a different choice of e) 

+ 2 H|Vu ||2 + v\\u\l <2 (||/|2+J7l 


where 


In particular, 


Vi ■= Vo + x 

2 La 


l«l- 


(14) 

(15) 


(16) 


(17) 


(18) 


^Ml + v\u\l < 2 (l/llz + »7i) ■ 

Using Gronwall’s inequality, we end up with the following uniform bound with 
respect to time 

Ht)\l< [nolle-’'*+ 1 (mlo,2+Vi)- (19) 

Now, integrating ED with respect to time, we obtain 

2V So H Vu ( s )ll2 ds + a J Q M s )ll2a+2(n) ds ^ ll^ola + 2 (ll/llL,2 +Vo) t. (20) 

Obviously, this inequality gives the required estimates for the existence of weak 
solutions for system 0,0. Furthermore, ED shows that 

limsup \u{t)\ 2 2 < - (||/||L ,2 + Vi) ■ 
t —>■+00 v \ ' 
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2.3. Gradient Estimate: Next, we show that Vu £ L 2 (0,T; L 2 (tt)) and that 
A u £ L 2 (0,T; For this purpose, we multiply the first equation of (HJ by 

—Am and integrate over H to get 

2^W S/U \\ 2 2+ 1 '\\ AU \\ 2 2+ (~ Au ) dx + a ( 1 + 2a ) II M“ Vu||2 

+b (1 + 2/3) | \uf Vu\l < i ||/||L,2 + \ II M\ 2 2 . (21) 

To obtain this inequality, we used the fact that f Q |m| 2 “m • (-Au)dx = (1 + 
2a) |||m|“ Vu|| 2 - Now, using Caucliy-Schwarz, Holder and Young inequalities and 
assuming a > 1 , we can write 


/ (it • V) it • {—Au)dx < / lullVitl^lVitl 1 “ | A u| dx 
J n J n 

< IIM |Vti|“| L 2 a (fl) || | V m| 1_ “ | 


(n) 


IIA i 


<i-||H“|VM||||||VM||f-^ + eo||AM||^ 


— T~ H u l a V u ll2 + ~T~~ 

4e 0 z 4e 0 


|V«|^+co \\^u\\l 


II2 T II ^ U II2 • 

In the sequel, we will need as well the following estimate 


\\\u\PVu\\ 2 2 <\\\u\ a Vu\\? |V«ir-" < e|H a Vu|® + 

Gathering both estimates, we get 


_ , f> 

1 




j t \S7u\l + u\Au\l+a{l + 2a)\\u\ a \S7u\\l < ^ I/I 2 + M ||V it^, 


where we set 

Ai := 2 


!/“a(l + 2a) 
.2(1 + l{t><o})_ 

In particular, we have 


4 /3 [H(l + 2/3)] c 
{b<0} [ [a(l + 2a)]/ 3 


l/la+^i IIVitla- 

Integrating (E5ll with respect to time, and using estimate (12U1) leads to 

|Vu(t)|2<I [(2 + Ai) 1 / 1 ^ 2 +^%] t +^|«ol 2 + |Vnolla- 
for all t £ [0,T]. Now, we integrate (l24l) with respect to time, and obtain 

v f | A zi(s)I 2 ds + a (1 + 2a) f |||it(s)|“ V ix(s)|| 2 ds 

J 0 J 0 


< 


1 r 


(2 + ^i) 


ll 00,2 


( 22 ) 

(23) 

(24) 


(25) 

(26) 


" 2 ’-Ai i?o t + ^IMa + |Vuo||a- 


(27) 


In summary, uni, and (1271) ') show that u £ T 2 oc (R + ; Fl 2 (fl)). In particular, we 
conclude that A u £ i 2 oc (R + ; i 2 (H)) and u £ L 2 oc (R + ; L°°(f2)), thanks to Sobolev 
inequality. Before we show the continuous dependence on the initial data, the 
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uniqueness, and the fact that dtu e ifoc(R + ;L 2 (ri)), we first improve the bound 

(EH). 

2.4. Uniform Gradient Estimate: Next, we show that Vu belongs to L°°(M + ; L 2 (£l)), 
that is ||Vu|| 2 is uniformly bounded with respect to time. More precisely, we claim 


Proposition 2.1. Let f £ L°°(]R + ; H) and ug £ V. Assume a > 1,0 < (3 < a,a > 
0, and b £ R. Then, the solutions of (CUD satisfy 


l|Vu(t)|| 2 < K := 


Ki if 0 < t < 1 , 

K 2 if t > 1 , 


where K i and K 2 are given in (ld2l) and (Iddl) . below, and A\ is defined in Theorem 

rn 


Proof. We start with the bound (l26l) . On the one hand, for all t £ [0,1], we have 

||Vu(t)||^ < i [(2 +Mr) l/lL.2+^%] +^M2 + |V«o|2- (28) 

On the other hand, for all t > 1, we consider s £ [t — l,t], and integrate (05l) over 
[s, t] to get 

\\Vu(t)\\l<\\Vu(s)\\l + - \\f\\l o 2 {t-s)+Ai I | V u(t) 1 2 dr 

J S 


V 


<lV<*)ll + l\fL,2 + Ai 


IIV u(t)\\ 2 dr. 


(29) 


't-1 


Notice that integrating (1171) over [t — l,f], and using (IT51l we obtain 


hoik 


2v It x II v u ( s )ll2< 2 (|/|L,2 + m) + IK*-Oil 


< 2 < ( 1 


III,2 + ^ + 77 




(30) 


Plugging (13U1) into (1551) . we get for all s £ [t — 1, t] such that t > 1, 


IIV u(t)\ 2 < |V w(s )|| 2 + —- e 


2 ^1 


1 

+ - 

V 


2 + A\ 


1 

1 + - 
v 


ho | \ 

l/lL.2+^1 


Vo H- 

V J 


Now, we integrate this inequality with respect to s over [t — 1 ,t\ and obtain 


II v u{t)\l < J* ^ IIV u(s )|| 2 ds + ^ e-^- 1 ) lluoll 


+ 


2 + Ai 


1 + 


ll/ll 


oo,2 


Ai 


Vo 


Vi 


<_ (l +A)e -Kt-i) | Uo | 


2 + (A± + 1 ) 


1 + 


II 00,2 


+ (Ai +1) 


Vo 


Vi 


(31) 


Eventually, (051) and m show that 

|V u(t)\ 2 < K\, for all 0 < t < 1, 
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with 


K i := 


A 


{2 + Ai) 11 / 1 ^ 0,2 + -^i Vo + ^ ll u o ||2 + ||Vuo| 


and 

with 

K 2 


\\Vu(t)\\ 2 <K 2 , for all t > 1, 


(32) 


2 + (A\ + 1) 


1 + 


loo,2 + ("4l + 1) 

1 

+ 2 ^ 


Vo 


Vi 


(1 + A\) IIMoII 2 


■ (33) 

□ 


It is worth noticing that this uniform bound is essential for establishing the existence 
of a nonempty compact global attractor of system ©,© when / is independent 
of time, and we refer to [20]. As we will see below, this bound will be crucial in 
showing the convergence of the solution of the data assimilation algorithm to the 
solution of (fTl2l) . Eventually, observe that using m, we obtain clearly the lim sup 
bound of Theorem O 

Now, integrating m, leads to the following inequality for alH > 0 


rt +1 


| Au(s )|| 2 ds + a (1 + 2 a) 


r t+i 


|u(s)r|Vu(s)|||^ 2 ds 


<-|/loo,2 + (^l+l)^ • 


(34) 


Now, we show that if the initial data Uq £ V fl L 2a+2 (fl), then we have u £ 
-k£c(® + ; L 2a+2 (Q)) and d t u £ L 2 oc (R + -,L 2 (fl)). For this purpose, we multiply the 
first equation of © by d t u and integrate over fl to get 


f \d t u\ 2 dx + f d t uu ■ X7udx + ^ f |Vit| 2 da; + 
Jn Jn 2 dt J n 


2 o :+2 


dx 


2a + 2 dt Jq 

f \u 2 P +2 dx = f fd t udx (35) 
Jn Jn 


2/3 + 2 dt Jq 

In particular, thanks to Cauchy-Schwarz and Young’s inequalities, we can write 


f \d t u\ 2 dx + f d t uu ■ Vudx + ^ [ |Vit | 2 dec - 

Jn Jn 2 dt J n 


2a+2 


2a + 2 dt 


dx 


2/3 + 2 


j t Jju 2 ? +2 dx< l/l 


(36) 
















3D BRINKMAN-FORCHHEIMER-EXTENDED DARCY MODEL 


13 


Furthermore, we have 

J ^ d t uu-Vudx< ^\\d t u\\ 2 2 + K^uf^ 

< i ||S Hl 2 +(^ + ^) \uf m +v\\Au\\l 

< ^;ll^t M ll2 + -^ ^ ^I12 + - (ll/llL,2 + ^l) + ^ ll^ U ll2 

: = \\\dtu\\l +v\Au\\l +P 2 (37) 

Thus, we can write now 


1 

4 


2 dt 


||Vu||( 


“K 2 dt 


Ml 


2 o ;+2 

2 q ;+2 


5 d 
2/3 + 2 dt 


2/3+2 

2/3+2 


< I/||Au||2 + |/||2 + P2- 


Integrating this inequality with respect to time and using (fT5l) (if b < 0) and (f27]l . 
we obtain clearly 

«e^ c (l + ;L 2 “ + 2 (ll)), $uGl£ c (R + ;I, 2 (n)), for all u 0 G V n L 2a+2 (n). 


Now we turn to the proof of the continuous dependence on initial data and the 
uniqueness. 

Remark 2.2. ^4s mentioned above, most of the calculation we presented are for¬ 
mal, but can be done rigorously by performing them at the Faedo-Galerkin level and 
passing to the limit. Nevertheless, the action of the equation on u makes sense, in 
particular since we have clearly d t u G L 2 oc (R + ; R _ 1 (n)) and u G L 2 (S. + ; H 1 (tt)) 
so that the function t e+ |u(t )|| 2 is absolutely continuos and the action ( dtu,u) = 
\j±\u{t)h.( see e.g. [301 EH]/- This operation is then rigorously defined when per¬ 
formed on the continuous equation. However, when we multiply by — A u, the action 
(dtu, — Au) is no longer equal to 1 4\ ||Vti(i)l 2 > except for initial data in L 2a+2 (fl) 
where we show that dtu and — Au are both in Tj 2 oc (R + ; L 2 (fl)). Therefore, the mul¬ 
tiplication by —Au and the integration over has to be done either on the Faedo- 
Galerkin level, or using a regularization procedure like convolution with a modifier, 
and pass to the limit. Eventually, when the initial data is in V f|I 2 “+ 2 ^), then 
8 t u G L 2 oc (R + ; L 2 (fl)) so that the multiplication of the equation by dtu makes sense. 
This remark holds true for all the systems we consider below. 

2.5. Continuous dependence on initial data and uniqueness of strong so¬ 
lutions. Let u and v be two solutions of system ([1]) , (El . corresponding to initial 
data uq and vq, respectively. Let w = u — v, then an easy calculation shows that w 
satisfies 


d t w — v A w + (w ■ V) it + (v ■ V) w + V (p u — p v ) + a (|u| 2q u — |u| 2 a u) 

+b (ju\ 2/3 u — \v\ 20 v ) = 0 . 
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The action of this equation on w leads to 

\j t \w{t)\l + Jj(w • V) u] ■ w + i/ | V HI 2 + a J n (I u | 2 “ m - \v\ 2a v) ■ w dx 

+b I (\u\ 2 ^u — \v\ 2 ^v) ■ w dx = 0. 

Jn 

It is well known (see, e.g ., [5]) that there exists a nonnegative constant kq = Kok) 
such that 


0 < no \u — v\ 2 (|u| + \v\f a < (\u\ 2a u — |c| 2q c) • {u — v). 
Moreover, we use the following classical inequality 

(\u\ 2 ?u - \v\ 2p v) < k 0 (|u| + H) 2 ' 3 kl- 
Now, using Holder, Young and Ladyzhenskaya inequalities, we can write 

J [(w -V)u]-w< | V u[ 2 \w\l< (a) 

< K il|Vii|| 2 ||w||| Hk 

< e I Vie 1 2 + ^ + pj- IVul^ Ma- 
Next, we use Holder and Young inequalities to get 


(\u\ 20 u-\v\ 20 v 


)-wdx<ko [ (kl + M) 2/3 |u>| 2 “ kl^ 1 “)da 

Jn 


x-P 


< K 0 | (|u| + |u|)“ k|| L ? (n) ||w|| L 2 ( “ a) 


Now, we use 


< eko | (kl + |w|)“ klk + k 0 

, optimize in e, and obtain 


^ , p 
1 


Ml 2 - 


— \\w(t)\\l + v I Vic 1 2 + (2 - l {b < 0 }) a Ko, II (kl + \v\) a kills 

(2k 0 |6|) c 


<L+^1k 4 + 1 


{b< 0 } 


with K as in Proposition 12.II In particular, it holds 


{an o y 


a — 0 


m;( 


dt 


\\w{ 


112 — 


< <u + ^K 4 +t 


{b<0} 


(2k 0 \b\Y 
{an o)P 


lk(<) lla- 


Using Gronwall’s inequality, we get 


H*)l2<H* = 0)l2e 

which obviously leads to the desired result. 


2 Jy v+1 ^ Ki+ 1{b<0} ] 3-5 } * 


(38) 

(39) 


(40) 


(41) 


2 

2 ’ 


Remark 2.3. In order to show the continuous dependence on the initial data and 
the uniqueness, we used the uniform bound given by Provosition \2.1\ Nevertheless, 
this fact can be shown without recourse to this bound by using the local integrability 
(in time) of\S/u(t)\^ given by (fl9l) and the estimate (E^l) . The proof we presented 
using the uniform bound on ||Vit(t) I 2 simpler. 
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Observe that ifuo,vo £ V nL 2 a+ 2 (ff), then u,v £ Lf oc (R + ; L 2 (fl)) 


Remark 2.4. 

so that dtw £ 

w makes sense rigorously. However, if the initial data uq and vq are only in V, 
then dtw £ L 2 nc (R + ; and we have to consider the action of dtw on w £ 


L 2 oc (R + ; L 2 (VL)) and the multiplication of the equation of dt:w by 


L 2 oc (R+-,H\n)). 


3. Proof of Theorem 11.41 

In this section, we prove Theorem 11.41 by seeking for a special solution of system 
Jl]),(0 °f the form u(x,y, z,t) = (4>(z,t), 0, 0), with </>(0, t ) = <f(L, t) = 0, and p = 0. 
Plugging these requirements in equation we find that 4>{z,t ) must satisfy 


d t </> — vd 2 <j> + a\(f\ 2a (j) + b\(/)\ 2 ^(j) = 0 , 

4{o,t) = <f>{L,t) = o. 


(42) 


For any smooth initial data, equation (1421 ) has short time existence and uniqueness 
of smooth enough solutions. In the sequel, we will show that for a special set of 
initial data <j)(z, 0) = 4>o{z), the solution of equation (l42l) blows up in finite time. 
For this purpose, we need the following 

Lemma 3.1. Let a, b £ R, 4>o(z) > 0, be smooth enough initial data , and 4>{z,t) be 
the corresponding solution of (l42l) . Let [0,T) be the maximal interval of existence 
for <p(z,t). Then, <f>(z,t) > 0 for all t £ [0,T). 


Proof. Let us denote by </> = max{0, —</)}. We will show that =0 for all 

t £ [0,T). We multiply (l42l) by </>“ and integrate over (0 ,L) and obtain 

7 --j- f \(f>~\ 2 dz + u f \d z (j)~ |“ dz + a f \(j)~\ 2a+2 dz + b f \(p~ \ 213+2 dz = 0. 

2 dt J Q J Q J Q J Q 

This equality implies that for all t £ [0,T) 
d 


dt 


I rt 


L 2 (0,L) ^ 2 


(H linii»(0,L) + |fe|«nif» ( 0,L)) in 


- II 2 

L 2 (0,L)~ 


Therefore, by Gronwall’s inequality, we have 


lr(t)\ L >(0,L) < e 2/ o( |a|W " (r)lli “(o,r) + H l^ ( r ) ir- ( o, t ,)^^-(0)|2 2(0ji)- 

Since \\tp~ (0) 1^2 (-q l ) = 0, and the exponential term is finite, we infer that \\4>~(t) 1 ^ 2(0 l) 
0 for all t £ [0,T). Consequently, <p(z,t) > 0 for all t £ [0,T). □ 


Now, let (fo (z) > 0 be smooth enough to be chosen later. Let T > 0 be the 
maximal time of existence of the solution </> of equation (|42[) . corresponding to the 
initial data (f> o. Thanks to Lemma lTll f>(z, t) > 0, for all t £ [0, T). We aim to show 
that T < +oo. We proceed by contradiction and assume that T = +oo. Then, 
we multiply equation (l42l) by sin ^ > 0, and integrate over (0, L), to obtain, after 
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integration by parts, 


d_ 

dt 


l‘ ^ 77 Z / 7T \ ^ f ^ 7T Z 

/ (j>(z, t) sin —— dz + v ( — ) / <f>(z, t) sin —— dz 

Jo L 'LJ J 0 L 


+ CL 



\4>(z, t) | 2a sin ^-^dz + b 

Lj 



\4>{z, t) | 2,9 sin ^ dz = 0. 

(43) 


Now, we need the following basic fact 

Lemma 3.2. Let 'ip £ L 1 (0,L) satisfying ip(z) > 0, a.e. on [0, L\; and let 7 > 1. 
Then, there exists a constant c 7 such that 



ip(z) sin —— dz ] < 0 ^ / %/>(z) 7 sin —— dz. 

L J Jo L 


Proof. The proof is a straightforward calculation. Indeed, we have 



The proof is achieved by setting c 7 


(fr 1 - 


□ 


Next, we introduce the notation m(t) := / Q L f>(z, t) sin dz > 0. Since (j>(z, t) > 0, 
by applying Lemma 15721 to equation (H51) and using the fact that a < 0, we obtain 

4- mi J) > ~ v (t) 2 m (J) + rn(t) 2a+1 -— m(t) 2/3+1 . 

dt \LJ C2a+1 C 2P+1 

Since a > (3 > 0, there exists to* > 0, large enough, such that for all to > to* we 
have 

- 1 / (4) 2 m(t) + -^-m(t) 2a+1 m (t) 2a+1 . 

'LJ C2ct+1 C2/3+1 4C2a+l 

Now, we choose to(0) > 2to*. Therefore, by the continuity of m(t), one has for 
short positive time that m(t) > to*. From the above, we have for short time 

^m(t)>Am(tr 1 . (44) 

dt 4c 2q +i 

Thus, we are able to conclude that rn(t) > m* for all t > 0 and as a result, 
(l44l) also holds for all t > 0. Since a > 0, then by integrating (l44l) we obtain 
m(t) > to(0) (1 — a 2 J°^ + i to(Q) 2q t)~^. Hence, there exists a finite time T* such 
that limsupj^j.* m(t) = + 00 . This shows the finite time blow-up and concludes 
the proof of Theorem 11.41 
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4. Proof of Theorems EH 

The proof follows the same line of the proof of Theorems 11.11 and 1 1.2 1 with the extra 
difficulty that consists in dealing with the difference I^(rt) From this point 

onward, we will no longer split the proofs depending on the sign of b and use the 
characteristic function H{b < o} f° r shortness. 


4.1. Local-in-time solutions. Let a be a strong solution of system (|T|) , ([? )) ■ our 
goal is to prove the existence and uniqueness of strong solutions to the following 
system 


9 f w-i/AB + (a-V)t) + V(; + a \v\ 2a v + b |u| 2/3 v = <p — 


Si : 


V ■ v = 0, v\ t =a = v 0 , 

v(x + L, y, z, t) = v(x, y + L,z,t) = v(x, y,z + L,t)= v(x, y, z, t), 


q{x + L, y, z, t ) = q(x, y + L,z,t)= p{x, y,z + L,t)= q{x, y, z , t), 


where we set ip := f + pih(u). Before going further, let us recall that u £ 
C°(R + ;V), by Theorem 11.21 so that thanks to the triangular inequality and ©, 
we have 


\ x li u ) 1 2 < \u-l\(u)\ 2 + \u\ 2 < y/co h ||Vtt| 2 + |«| 2 . (45) 

In particular, this shows that I^(u) £ C°(K. + ;H), which in turn implies that ip £ 
C'°(R + ; H), and that for all t > 0 it holds thanks to (TT9l) and Proposition 12.II 

IMI 2 < := l/l^ + py/c^hK + p ^||u 0 | 2 + ^ (||/||L , 2 + 7 ? 1 )) • ( 46 ) 


As in section EH we start by showing the existence of local (in time) solution 
to the system <Si by the mean of Faedo-Galerkin approximation procedure. More 
precisely, as in m, we expand v as follows 

m 

v m (x, t)=Y^ 9 l m(t ) Wi(ar) £ V m , (47) 

k=1 


9 m /n 


Wi dx , for i = 1 , 2 , ..., to , are the unknown coefficients, solving the 


following system of ordinary differential equations in V 
d 


dt 


Vm-Wkdx + v / Vv m • VW fc dx+ (v m • V)v m • Wfc dx 

J £2 J Q, J £2 

W k dx + b |v m | 2/3 v m • y\4 dx 
J n 


|v m | 2 “v m 


= / (p ■ Wfe dx — y / X£(v m ) • Wkdx, 
Jn Jn 

9m(t = 0) = /vo • W k dx, 


(48) 

(49) 


for all k = l,... ,m. The vector field in system (TT71) . (TTfTl) being obviously locally 
Lipschitz, the system admits a unique solution g^t) £ C 1 ([0,T m ]), for some small 
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interval [0, T m ] C [0, T], As in the previous section, we will perform formal estimates 
using system <Si instead of the Fadeo-Galerkin system for simplicity. 

4.2. Global-in-time solutions. We multiply the first equation of system <Si by v 
and integrate over Q. Using Cauchy-Schwarz and Young inequalities we obtain 

t; 4 Ml + v l Vu l2 + a Mla+l + b Mlp+l 


2 dt 


< - Ml - 7 M Ml + M / (v(t, x) - iMt, x))) ■ v(t, x ) dx 
t L 4 Jn 


< Ml - \vMl + mc 0 /i 2 IVc 12 . 

Now, thanks to msD, (m, and the fact that 2/j.cq h 2 < v 1 we obtain 


IMI 2 + M Ik II2 + ^ II V4 ’ll 2 + (2 - 1 {6<0}) a Mla+l < 2 


dt 


M 


?72 


(50) 


where we set 772 := 1{6<o} H 


2|b| 


3+1 


|f2|. In particular, we have 


M 


— IIW1 2 + At Ik II 2 < 2 ( — + 772 




and by the virtue of Gronwall’s inequality we get 


lk(t)||^<lKll2e-^ + - (- + %)• 


M 


1,2 ^ „2 t + 2 

Ai V M 

Notice that integrating (f50j) with respect to time, we obtain 

r t r t 


(51) 


P / ||v v(s )||2 ds + (2 — l {b<0} ) a f |k(s)l 2 a +2 ds 
J 0 Jo 


< Mil 


M 


+ m I t. 


(52) 


Now, we find a bound on ||Vu| 2 . For this purpose, we multiply the first equation 
of system <Si by —Av and integrate over 12 (see Remark 12.2D . We obtain 


l^j:\\Vv\\l + i'\\Av\\l + J^(v- V)v ■ (-Av)dx + a(l + 2a)\\\v\ a Vv\\l 

+ b (1 + 2/3) Ikk Vv\\l = f ip ■ (-Av)dx-fj, f X^v) ■ (—A v)dx. 

J Si J SI 

Now, using Cauchy-Schwarz inequality and dSD , we can write 

[ lk(v) ■ {~Av)dx = n [ (v-Xl(v)) ■ {-Av) dx ~n [ (-A«) 

J SI J £1 J SI 


< 


At 2 co h 2 

ir - 


- ti) | V v\\ 2 + e | Av|| 2 • 


dx 


(53) 


Using this, and the fact that 2/7 co h 2 < v, and proceeding as in the previous section, 
by using (l22l) and (l23l) . with different optimization in the e's, we get 

4 I1 2 + * |At,|* + a (1 + 2a) \\v\ a V«| \ \ Ml + A 2 |V«|*, 


(54) 
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where 

= h f [ v a a (1 + 2a) 
2 + \[3“(l + l {b<0} ) 

In particular, we have 


46<0} 


4^[|6|(1 + 2/3)] c 


[a(l + 2cl)]P 


d „ — 1,2 , 2 n n 2 


-||V<<-|M|^ + .4 2 1Vv i 2 - 

Integrating this inequality with respect to time, and using (1521) we obtain 


A 


II Vw(i )||2 < — \vo\\; + |Vu 0 1 2 + - M 1 + ^ + A 2 V 2 )t. 


Ad 


Moreover, integrating (15H) and using (15^1) . we get 


(55) 


(56) 


" / \\Av(s)f 2 ds + a(l+2a) / || v ( a )|“ V«(a)|^ 


A: 


— — II ^ 0 II 2 + II ^ v 0 1 2 + 7 ( M ( 1 + — ) + -4-2 »72 ) t. 


Ad 


The global existence of solutions to <Si follows as for Theorem 11.21 Let us mention 
that the calculation we performed in this proof is formal. Nevertheless, it can be 
done rigorously using the Faedo-Galerkin system, especially when multiplying the 
first equation of system 5i by — A v and integrating over Q and we refer the reader 
to Remark |2.21 


Before going further, we show the following 


Proposition 4.1. Let f £ L°°(R + ;H), u the strong solution of system 0.0 
given by Theorem, \1.2\ vo £ V, a > 1, 0 < /3 < a, a > 0 and b £ K. If 2 ^ co h 2 < fi, 
then, there exists P > 0 such that the solutions of m satisfy 


|Vu(t )| 2 < P := 


Pi if 0 < t < 1 , 
P 2 if t > 1 , 


where 


Pi := 


Pd ~ 


— IMI 2 + ||Vuo || 2 d— -) + v4 2 ? 7 2 

v v \ \ /i ' 


1 


- (A 2 + 1) | uq 1 2 + - < M + (A 2 + 1) 1 + --b in 


1 


T 


M 


T 


and A 2 is defined in (1551) and g 2 := l{ b <o} 216| 


/3 + 1 



a 


N- 


Proof. The proof is in the same spirit as the one of Proposition 12.11 On the one 
side, if t £ [0,1], then we use (15B1) to get 

||Vv(f )||2 < — ||uo ||2 + IIVU 0 I 2 + - ^1 + —^ + A 2 t? 2 ^ . 

On the opposite side, let t > 1 and introduce s £ [t — l,t]. Using (l5il) . we have 

o M ft 

IIV v(t)\\l < |Vu(s )||2 +- + A 2 / |V«(r)|“dr. (57) 

v Jt-i 
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Now, we integrate (l50l) over [t — 1 ,t] and use (1511) to get 

v £ | V v(s)||2 rfs < 2 ^1 + ^+e _M(t_1) l^oll 


Therefore, © becomes now 


\\S7 v (t)\\ 2 2 <\\\7v(s)W 2 2 + ^e-^ 


+ -{M + A 2 1 + 


M 


- V2 


Next, we integrate this inequality with respect to s to obtain 


ivw(*)ii < r ivu( S )i^ s 

Jt-i 


— e-^\vo\l 


M 


4— < M + A 2 ( 1 H— [-+V2 






< l (A 2 + 1) e-^-V |«o|= + 1 |m + (A 2 + 1) (f + 


M 


- rj2 


This leads to the desired result and achieves the proof of Proposition l4.il □ 

Now, we show the fact that v G I/^ C (M + ;L 2 a+ 2 (fl)) and dtv G L 2 oc (M + ,L 2 (fl)) 
if vo G V fl L 2 “ + 2 (fi). Actually, it follows the same argument as in the proof of 
Theorem 11.21 in section [2] and the extra term to deal with is the following 

J^ll(v)d t vdx < e\d t v\ 2 2 + \li{v)\\. 

Now, optimizing in e, the term e || 9* v || 2 can be absorbed in the left hand side of the 
equivalent expression of (IQ) for v. Eventually, inequality (Pf5l) shows clearly that 
P »|| 2 is integrable in time which leads to the desired result. 


4.3. Continuous dependence on the initial data and uniqueness of solu¬ 
tions. Equivalently to (JT] [ 2 ]), before we show the continuous dependence on the 
initial data and the uniqueness of such solutions, let us prove the following Now, 
we turn to the proof of the continuous dependence of the solutions v(t), obtained 
in the previous section, on the initial data and their uniqueness. We mention that 
Remark no is still valid in the case of system S\. Let a be a strong solution of 
system 0,(11 and Vi and v 2 two solutions of 00 . Furthermore, let w = v\ — v 2 , 
then w satisfies 

dt w — v A w + (w ■ V) vi + {v 2 • V) w + V (p Vl — p V2 ) + a (|vi| 2 %i — \v 2 \ 2a v 2 ) 

+b (|vi| 2/3 ^i - \v 2 \ 2/3 v 2 ) = -pl l h {w). 

(58) 

Now, clearly w G L 2 (R + ; fJ 1 (fl)), and since d t V\,d t v 2 G Lf oc (R + ; then 

we have also d t w G Lf oc (R + ; H~ 1 ( fl)). Therefore, the action of the equation on w 
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leads to 

1 d . 

— — \\ w 

2 dt 


We will need 


+ f [(w • V) Vi] ■ wdx + v |V tofl 2 + a f (\v\\ 2a vi — \v 2 ,\ 2a V 2 ) ■ wdx 
J n J n 

+b / (}vi\ 2l3 vi — \'V 2 \ 2!3 V 2 )-wdx = —fi / l\{w)-wdx. 

J n J o 

(59) 

J ^ X h( w ) ■ wdx < tlQlL |V w\\l + (fie - fi) \wf 2 , 

, and optimizing in the e's, we obtain 


w(t) \\l, (60) 


Now, gathering this inequality and 

^ \\w(t)f 2 + v I' Vwf 2 + (2 - 1 {6<0 }) a k 0 | (|u| + |„|)“ M I 2 


A-»+\ v - 


2 7 k? 


P 4 


1 


[(2k 0 |6|)“1^ 
{b< ° } [ (a« 0 )^ . 


Gronwall’s inequality leads to the desired result. In particular, we have 

r t+T r t+T 


| A v(s )|| 2 ds + a (1 + 2a) 


Ns)r V u(s )|| i2 ds 


<-|Mll.2 + (-42 + l)P 2 . 


(61) 


Observe that if no £ V fl L 2a+2 (tt), then d t Vi,d t V 2 £ Lf oc (R + ;L 2 (H)) so that 
dtw £ Lf oc (R + ;L 2 (Q)) and the inner L 2 inner product of dtw with w is well defined 
and no regularization is needed. 

5. Proof of Theorem 11.61 


5.1. Convergence in L 2 norm. In this section, we show the convergence of the 
difference u(t) — v{t) to zero as t —> +00 in the L 2 norm where u(t) is a strong 
solution of system 0.0 ensured by Theorems II. 1 1 and 11 .2 1 and v(t) the solution of 
system 5i with an interpolant I\. This will be achieved by suitable assumptions 
and the smallness of h and that fi is large enough. Let w = it — v, then subtracting 
the first equation of system <Si from the first equation of system 0.0 leads to 
equation (1551) (with iq replaced by u and v 2 replaced by v). We multiply this 
equation by w and integrate over 12 to get 

||iix(t)I2 + v ||Vu>||2 + [ [(w • V) u] ■ wdx + a f (\u\ 2a u — \v\ 2a v) ■ w dx 
J n J n 


2 dt 
+ b 


(I 


it| 2/3 w - |d 2/3 u 


) • wdx = —fi / Zh(w) ■ w dx. 

J n 


Therefore, (15U1) still holds true with P replaced by K. Observe that I\ is indepen¬ 
dent of fi. Therefore if we choose fi large enough such that 


1 o7*.8 


f ‘ >2 2 F+ „ 


2'kX , 
i- TS 4 


r(2«o|6|)»]' 

K +1{b<0} [ (a Ko )? \ 


then, using Gronwall’s inequality we obtain 

\\w(t)l<\\w(t = 0)\\ 2 2 e-^. 
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This shows the convergence of u(t) — v(t) to zero, as t —> +oo, in the L 2 norm. 


5.2. Convergence in H 1 norm. In this section, we show the last assertion of 
Theorem H31 namely the convergence of the difference u(t ) — v(t) to zero, as t —> 
+oo, in the H 1 norm where u{t ) is a strong solution of system dT|) , (j2j ensured by 
Theorems 11.11 and 1 1. 21 and v(t) the solution of system S±. Now, since we restrict the 
range of a to 1 < a < 2, therefore we have d t u,d t v £ Lf oc (R + ; L 2 (tt)), therefore 
d t w £ L 2 oc (R + ; L 2 (£l)) and —Aw £ Lf oc (M + , L 2 (fl)) so that the L 2 inner product 
of (l58l) with — Aw leads to 


| V I 2 + v | A w \ 2 2 + f [(«; • V) u] • (—A w) dx + f [(u ■ V) w] ■ (—A w) dx 
2 at ~ Jq J 0 

+ a f (\u\ 2a u — \v\ 2a v) ■ (—Aw) dx + b j (\u\ 2 ^u — \v\ 2 ^v) ■ (—Aw) dx 

J n J n 

= — /j j X^w) ■ (—A w) dx. (62) 

Ja 

We need the following estimate 


[O' V)ti] • (-A w)dx < |Vu|| 2 ||«;| 0O ||A«;|| 2 < n 2 ||Vu|| 2 Ml Id \w\ 2 H2 ||Aw|| 


< k 2 |V«| 2 |w||ri \\w\\^ 2 < 


-4 IIV it 1 2 


Mm + e |Awt| 2 . (63) 


Also, we need 


[(u • V)w] ■ (-Aw)dx < ||Vw|| 2 | Aty|| 2 


< «2 ( ~ e M ff i + e \v\\ H2 ) |Vit?j| 2 |Auj|| 2 


< 


k\ 


1 


4ei V 16e 2 


Ibllffi |Viu| 2 + 2ei||A«;|| 


„ 2,2 

2 K 2 e 


4ei 


II At'1 2 |Vw|| 2 . (64) 


For later use (see below), it is crucial to avoid having the term ||Au|| 2 in the right 
hand side of this inequality. For this purpose, we use the triangular inequality to 
get 


||A-u1 2 ||Vw| 2 < 2 (||Aw1 2 + |AitI 2 ) |V-w 1 


< 2||Au|| 2 ||Vu)|| 2 + 4 
Therefore, inequality (l64l) can be replaced by 

f [(v • V) iy] • (—A w)dx < ( 2ei + 


II 1 2 ) ||Aw|| 2 . 


ei 


4ei V 16e 2 


|V«|1 + ||Vu|| 2 


. 2^2 


I Aw 1 2 


(65) 


Hffi|VH2 + v-|A«l2|VHS- (66) 


The difficulty in showing the convergence in the H 1 norm, compared to the L 2 (fl) 
one, consists in dealing with the power terms. Indeed, on the one side, there exists 











3D BRINKMAN-FORCHHEIMER-EXTENDED DARCY MODEL 


23 


Ki(a) such that 

f (\u\ 2 a u-\v\ 2 a v)-(-Aw)dx<K 4 (a) [ (|u| 2a + |i;| 2a ) |w|| — Aio| t2®, 

Jn Jn 

For simplicity, we set K 4 := Ki{a) > « 4 (/?). In particular, using this and Young’s 
inequality, we obtain 


(I 


u\ u — \v[ 


y) ■ (—A w)dx + b f (\u\ 2a u — \v\ 2 a v) ■ (— Aw) dx 
Jn 


<K 4 (a + e\b\) / (ju\ 2a + \v\ 2a ) \w\\ — Aw\dx + 2 n 4 \b\ 


^ , f> 
1 ' 


IIVid1 2 - (67) 


In order to deal with the right hand side of this inequality, we will need to restrict 
the range of a. In fact, on the one hand we have 


K 4 [ M 2 “|w|| — Au)| dx < K 4 |u|| 2 “ h||oo ||Azd ||2 
Jn 


< K4K 2 ||w||4“ Hjfi \\w\\^ 2 


< 


< 


< 


4 4 

^ l< ■+ 
h 


innuft+e 


WMm + e ll Aw lls 


1180 ; 


\\w\\ 2 Hl +e\\Aw\\l 


4 e 3 


Ikllffi+e||Au>|| 2 . ( 68 ) 


In order to pass from the third to the fourth line above, we used Holder inequality, 
which requires a < §. On the other hand, let a > §, then we can interpolate 4a 


between 6 and + 00 , that is 


K 4 


|u>|| — Aw| dx < 


A a 


\H\Z iVtula + e|Ato|j 


A .A 


K 2 t ^4 I. || 12 || |i4 (2a—3) iir, m 2 m a | 

< II“lie IIw||°o l|Vzo| 2 + e||Au;|| 


J 8 (a- 1 )) 12 4 


< 


4e J 


^ mT +3> mT-” iv^il+ e |4»l 


Now, we use the fact that 


M$ a+3) mi! 2q - 3) < 


IIH 2 


2a—3 
4 —2a 


hi 


2a + 3 
2 —a 

H 1 


+ e|M| 


H 2 i 


which makes sense only when 2{2a — 3) < 2, that is a < 2 to get 


k 4 / \u\ 
Jn 


|w|| — Atu| dx < 


_8(a-l) 12 A 


4 e 3 


2a — 3 
4 —2a 


2o:+3 

2 - 


I “Iff!" +SH 


H 1 


|Vt 


_8(a-l) 12 ,A 


I Azt || 2 || Vzd|| 2 


^IIAHIr 


(69) 


4e 3 II--112 II • “112 ' -II-112- 

Obviously estimates (l68l) and (l69l) hold if u is replaced by v and K replaced by 
P. However, since P depends on fj, (see Proposition 14.11 and (H51) b we will have to 
solve when 1 < a < | an inequality of type n > Const/4 Sa + Const 2 ^ 2 ( 2 “+ 3 ) + ... 
which does not have a solution most of the time. At least, in our case we cannot 
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exhibit precisely the solution’s range if it exists. The hypothesis on the initial data 
luolffi < A', ||t)o|| ff i < K is introduced in order to get over this difficulty. Indeed, 
by continuity of ||v(t)||#i, there exists a short time interval [0, T) such that for all 
t £ [0,T), it holds \v(t)y H1 < 3A' 2 . In the sequel, arguing by contradiction, we 
will show that, actually, we have T = +oo. Obviously, this assumption can be 
removed by assuming |Vuo|| 2 < K since ||u(t)|| 2 is bounded. We work with the first 
assumption for lightness of notation. 

In the sequel, we will need the following version of Gronwall’s Lemma for which we 
refer to [T9j 

Lemma 5.1. Let £(t) be an absolutely continuous and locally integrable function 
satisfying 4- £(i) + £(t) ((t) < 0. Assume that for a fixed s > 0 


pt+s pt+s 

liminf / £(s) ds > 5 and limsup / £~(s)ds 

t—>oo J t t—¥o o Jt 


< 00 , 


where S > 0 and f = max{—a,0}. Then, £(i) -A 0 exponentially, as t —> oo. 

Now, we assume that T is the the maximal finite time such that |c(t)||^i < 3A' 2 is 
satisfied. We separate the proof to two parts depending on the range of a. 

The case 1 < a < §: On the one hand, gathering the estimates used in section l4~2l 
with (I7ITTI RT8l) and (138114T1) . we obtain for all t £ [0,T) 

1 d 


— — |Vw(f)|| 2 < — v + £o + £3 + 2ei + 4ae 4 + 


p v 
2^ 


T 


3k| 

4ei 


2 4 e 2 


£2 


K 2 + 2Ki\b\ 


£2 2 r 
£i 


3/7 2 + K 


II Aw (i) 1 5 


rv 2 c 2 

^7 


|A«Wl2 |Vw(t)||^ 


4^ + eo + 4a e4 + ^l|nr 

e 3 


2e| 


340 :^- 80 : 


I< o 


\\w(t)\\ 


H D 


where, thanks to (fl9l) and Proposition 12. 11 we set 

(ii/iil,2 


im<< K 2 


M2+ 2 

V 


■m 


:= K. 


(70) 


Observe that K does not depend on fi. On the other hand, following paragraph 
ED we have 

1 d „ ,. s „2 are 0 


2 dt 


Mt)\i 


2-1 


{b< 0} 


+ M) Q w|| 2 < e 5 - v + 


8ee) 


|Vw(t)||( 


8 7^4 


+ | £5 + — M + 


k\K 


|6|k 0 


2|6|k 0 


ati 0 


IK*) II2- 


(71) 


Summing the last two inequalities, using the fact that 2p,coh 2 < u, and optimizing 
in the e's, we obtain for all t £ [0, T) 

— K*)ll jU + v IIAwfa + (2 - 1{(,<o}) a,n 0 ||(|u| + |v|)“ wf 2 


dt 


< 


S 3 — p, + max < <5i, 82 + 


2 5 (3K 2 + K 2 ) 


\Ml 




(72) 
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where 


. 3 2 7 k\ AT 4 


(2|6|k 0 ) c 


^2 = 


3 vK 2 


(aK 0 )P 
2 10 3 k%(K 2 +K 2 )K 2 


S 3 = 


2 6 (A ' 2 + K 2 ) 
2 7 k$K 4 5 


+ 4k 4 |&| 


915..4 4 8a 4 

Li Li A. 2 '*•'3 


2 4 i/3 

Now, we apply Lemma T5. II We write (1721) as follows 


\n\ 8a (s ia K Sa + Kl a ^j . 


+£WIMIffi < 0, 


with 


Thanks to (l34l) . it holds for all T > 0, 


£(s) -.= 11—63 — max ^ (Si, J 2 + 


2 5 (3 K 2 + K 2 ) 


|A«(*)|; 


V 

T 


rt+T 


\Au(s)\\lds<-\ff 2 + (A 1 + l)K 2 . 


Setting T = 1 leads to 


f t+1 - \ 

/ £(s) ds > /i — (S 3 — max < <Si, 62 + 


II 00,2 


+ v(A\ + 1) A 2 


2 5 v(3K 2 + K 2 ) 


Therefore, if we assume 


H > 263 + 2 max < (Si, <S 2 + 


2 ii/iiL,2 + v {Ai +1) a 2 


2 5 v(3K 2 + K 2 ) 


then, we have 


liminf [ £(s) ds > ^ > 0 and [ £(s) ds < — < + 00 . (75) 

t->oo J t 2 J t 2 

Thanks to Lemma 15.11 there exists a nonnegative constant 77 such that for all 
t e [0,T), it holds 


ft+i 


(73) 


(74) 


H\ 2 Hl <\\w(t = 0 )\\ 2 Hl e~ vt . 


(76) 


In particular, we infer that ||u;(t)|| ff i < 2K for all t £ [0,T). This implies that 
\\v(t )\\ H1 < 3 AT. We reach a contradiction with the fact that T is the maximal time 
for which ||u(f)||^i < 3 K 2 and shows that T = + 00 . Clearly, the estimate (l76l) 
shows the exponential convergence of v(t) toward u(t), as t —> + 00 . 


The case | < a < 2: The proof is similar to the previous one except the fact that 
now we use the inequality (1551) instead of (1551) . However, thanks (IClT) . we have for 
all T > 0, 

c\ 

\Av( s )\\lds<-Ml >2 + (A 2 + l)P 2 , 

The bound involves /r (through P) and this will lead to the same kind of problem 
we pointed out in the previous section about the powers of fx. To get over this 
difficulty, we will use different estimate than (1691) for v. Indeed, the only term one 
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has to handle is the first one of the the last line of inequality (l69l) . Using (l65l) . it is 
rather easy to see that 


ft 4 


J |u| 2 “|u;|| — Aw\ dx < 


K 


8(a-l) ..12 ,.4 

fv 3 fv 4 


4e 3 


,8(q:—1) 12 4 , 

>2 ^3 ^4 fc 0 


2 a-3 

4—2a 2a+3 


|V«|1 + e \Aw\l 


H\m a +*+2eo IIV w || 2 - (77) 


Let 6 := ^ k| 2 k\. Therefore, collecting all the estimates and using ( 1771 ) . we 

obtain 


1 d 2 aK 0 


{h<0} 


M + H)“H2< 


+ 4 e 5 + £0 + /4£6 — /Li + 


icfA' 4 K. 


"A' 4 + |&|« 0 


2 | 6 |« 0 ' 


a/^o 


/«2 e 2 2a0£ 7 \ 


^2 c 2 

£1 


,3 

e 4 


3/4 2 + A' 2 


M*)ll2- 


| Aii; (t) | r 


+ 


u ^ 3 ko 


2 4 e| 


£2 


A' 2 + 2k 4 |6| 


ad 

24 


3? 

^7 


2a — 3 
4 — 2a 


AT— + 3e 7 /v 2 


K 2- “ + A 


17 k 2 T-i 

65 ~ v + o7" + + £ 0 + 

5£6 e 3 


r K 2 e 2 a0(l + 2e 7 ) 


2e, 


24 


|V«j(f)|2 
|Vw(t )| 2 

||Au|rt |Vtxj(t)I 2 . 


Summing-up this inequality with (1711) and optimizing in the e's, we end up with 

d 
dt 


+ 4l A 4l2 + (2- l {h <0 })aKo |(H + |u|) a ic1 2 


< 


-fi + max jji, 82 + <5 3 lAulljjj ||ic| 


2 

H P 


where 


r . 2 (k? + 2 6 k 4 )A' 4 

<5i — 2zz -|---1- 2|6|/«o 


2|6|kq 


an 0 


= 


2 ■ 3 7 4 (3AT 2 + A' 2 )A' 2 vK 2 


2 6 3 6 0a 4 a0 


2 6 3 6 9a 4 a9 


3K 2 + K 2 

2535 a 2 05 ( 3/? 2 + / 4 2 )5 


+ AT 2 + 2k 4 |6| 


(3 


2q —3 
4 —2a 


~ 2q+3 

AT— 


a 2 05 (3 A' 2 + AT 2 ) 5 


- K 2 


K 2- “ + K 


+ p + 


2 7 4A ' 4 


<53 = 


2 6 3 6 (3A' 2 + AT 2 ) 


Ko + 0a 4 ( 1 + 


a 2 6^(3K 2 + K 2 )§ 
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As done previously, we apply Lemma T5. II bv setting 

£(s) := n — max |d 1; S 2 + h || Au|| 2 |. 

Now, we set T = 1 in (ITdl) and obtain 

£ £(s)ds > /x-maxjiJi, S 2 + ^ ll/llL,2 + C^i + 1 ) K 2 ^ 

In particular, the assumption on the size of /i, 

II > 2max /<5i, S 2 + S 3 (- \\f\\ 2 2 + (Ai + 1 )K‘ 


guarantees ca> and 


^ \\w(t)\\ 2 Hl +£(t) \\wf Hl < 0. 


(78) 


Eventually, Lemma 15.11 allows to reach the contradiction with the hypothesis that 
T < +oo exactly as done for the case 1 < a < §. As a matter of fact, we conclude 
to the exponential convergence of v(t) toward u(t), as t —> +oo, with respect to the 
H 1 norm. 


6. Proof of Theorem o 


6.1. Existence of solutions and convergence. In this section, we prove Theo¬ 
rem o We are concerned with the following system 

d t v — v Av + (v ■ V) v + V q + a \v\ 2a v + b |u| 2/3 v = f + p(X 2 (u) — X 2 (u)), 


S 2 ■ < 


V • v = 0, v\ t =o = vo, 

v(x + L, y, z, t ) = v(x, y + L,z,t) = v(x, y,z + L,t)= v(x, y, z, t ), 


q{x + L, y, z, t) = q{x , y + L,z,t)= p{x, y,z + L,t)= q{x, y, z, t). 

Instead of S 2 , let us denote by v = w + u and consider the following equivalent 
system for w 


dt w — v A w + (w • V) u + (v ■ V) w + V (p u — p v ) 


S„ 


+a (|u| 2 “u — |u| 2 “u) +b (\u\ 2 ^u — \v\ 2 ^v) = —pX%(w), 


V ■ w = 0, ui| t=0 =u 0 - vq, 

with the associated periodic boundary conditions. The local (in time) existence is 
readily obtained using Fadeo-Galerkin approximation (H51H51) with interpolant X 2 
instead of X\, and observing that in the left hand side v can be replaced by u — w 
and we leave the proof for the reader. Now, we establish the necessary a priori 
estimates for the global existence of solution to system S w which in turn implies 
the global existence of solutions of 1 S 2 . Again, we perform formal calculation using 
S w for simplicity. 
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In the sequel, we will use the following estimate which is the analogue of (1531) for 
I 2 instead of X 1 


-M / Zh(<P)- (-A tp)dx. = n / {-A<p)dx-n / <p- (-A ip)dx 

JQ JQ J 


< 


4 Co 


h 2 


4e 


— M ) |V^ 


/J, 2 Ci ft 4 

iT - 


|Ap| 


2 > 


(79) 


for all G H 2 (fl). In Theorem 11.71 we assume that ||uo||^i < A and ||uo||//i < K so 
that |iu(i = 0)| H i < 2 K. Again, by continuity of ||w(t)|| fl -i, there exists a short time 
interval [0,T) such that for all t G [0,T), it holds |to(i)||#i < 4A 2 . In particular, 
it holds |u(t)| ff i < 3A' for all t G [0, T). In the sequel, arguing by contradiction, 
we will show that, actually, we have T = +oo. Therefore, we assume that T is the 
the maximal finite time such that \\v(t)f m < 9I \ 2 is satisfied. On the one hand, 
we multiply the first equation of system S w by w and integrate over f2. Thanks to 
459} with I 2 instead of , u instead of iq, v instead of , and Hi, we get 




an o 

~T~ 


, | \a. 1,2 , 4o4 

M ) w h < 


4ei 


| Aw || 2 + 


p 2 c 0 h 2 

£o 4 -:- v 


4ei 


£o + £i H—y | Vit1 2 + 1{6<o} 


[ 2^(|6|A 0 )“1 

1 -i 

**-P 

(awo)' 3 

~ V 




|VH | 2 

(80) 


On the other hand, we use (l63l [Ml , [67} [68} and (l79l) and get for all 1 < a < § 


14- IIVHI 2 < 

2 dt " 112 “ 


4 a : 4 


e2 


+ 


+ 


/i 2 coh 2 

4e 5 


2«e 6 e 8 { a 44 “ K 4 £ 6 , 0 |g 

+ A 3 ^ 

e 7 4e 7 eg 

9 k\K 2 ( 1 


(A 8 “ + (3A) 8a ) 


HI 


m 


4C3 


V 16e? 


■ 2 r 4 | 6 | 


2 '2 


tine 


2 C 4 

2 e 3 


II Aw II2 


—V + £2 + 2e 3 + £5 + 2a«4fi6 T 


2(16668 + p 2 Ci 


h 4 


£7 


465 


^-2 C 4 

£3 


( K 2 + 9A 2 ) 


IIVzdI^ 

||Aw;|| 2 . 

(81) 


Now, using 4EZ}, the fact that 15/rmax{co, yTT} h 2 < u, and optimizing in the e's, 
we obtain 


^\\w(t )\\ 2 H1 +v\\Aw(t )\ 2 2 


< 


S 3 — P + max < 62 , S 3 + 


15 (A 2 + 9K 2 ) 


II Aw |: 


4)1 


2 


(82) 
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where 


4 2 • 15 3 fto/v 4 

<5 3 := —m - 1 


2 3 15 


c „ 2 • 15 3 k?A' 4 ^ 

So '■— 2i/ H-r-h 1 


<5i := 4 k4|6| 


\b\ 


v u 

/3 


{b<0} 


(2|6|« 0 ) Q 

(an oY 


9-1 5 2 4(K 2 +9K 2 )K 2 
2 5 is 3 


vK 2 


10 (K 2 + 9 K 2 ) 


Next, we apply Lemma lS.ll bv setting 

€(s) := fi — S 3 — max |j 2 , S 3 + 
Thanks to (l34l) . setting T = 1 leads to 


II An 1S 


ri+l r 

J £(s) ds > n — S 3 — max < S 2l S 3 + 


15 (Jv 2 + 9 K 2 ) 

f +u(Al + 1 ) K 2 ' 


This, if we assume 


H > 2 S 3 + 2 max < S 2 , S 3 + 


15(A' 2 + 9A' 2 ) 


+ v(jA 1 + 1 ) K 2 


15(A" 2 + 9A" 2 ) 


(83) 


then (1751) holds true. In particular, thanks to Lemma [5Tl the contradiction leading 
to the fact that T = +00 is obtained as in section 15.21 Thus, we conclude to 
the global existence of solution to S w which in turn infers the global existence of 
solution to S 2 since v = u — w. Moreover, the convergence at exponential rate of 
v(t) toward u(t), as t goes to + 00 , in the H 1 norm follows. 

Now, we handle the case | < a < 2. Using (J631 [M], EElI [Ml 113 and (1791) . we get 
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(84) 
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As in the previous case, using ([87]) . and the fact that 15/xmax{co, y/ci} h 2 < v, and 
optimizing in the e's, we obtain 






< 


where 


d 3 - i-i + max \ S 2 , S 3 + — ( 1 + 


2 (K 2 + 3 AT 2 ) 


WMt 


83 := — v + 


2 • 15 3 4AT 4 


15 


82 := 2 v 


+ 


A3 K 8a 


v J 
3^8 


4 6 15 12 K2^ 8Q 7 ^K3 8 k\ 2 


2•15Wf K 


V u 


M*)l 


H 1 ’ 


:= 4«4|6| 


, \( 2 \b\a 0 r]^ 

+ i{b<0} ’ 

6 n\{K 2 + 9K 2 )K 2 vK 2 


6 • 5 3 (AT 2 + 9 K 2 ) 


(85) 


The rest of the proof is the same as in the previous case and leads to the following 
condition on /x 


/x > 2^3 + 2 max 



r v 
<53+ 15 



11/111,2 + v iAi + 1) A 
2v{K 2 + 9.A 2 ) 



( 86 ) 


In the previous section, we used the fact that Vrt, Vw £ A°°(R + ;A 2 ) in order to 
show that dtu,dtv £ A 2 oc (R + ; L 2 (fl)) if the initial data is in L 2 a+ 2 (Vl). In this 
section, we need to use a different estimate instead of (l37l) . Indeed, thanks to (l22l) , 
one can see that the following inequality holds 


dtv v ■ Vx> dx 


< 


4eo 


! Vu \l 


4eo 


||Vu||2 + eo \dtv\l- 


Clearly, the term e 0 can be absorbed in the left hand side of the equivalent 

expression of (1361) for v and the remaining terms are obviously integrable in time. 
The other new term compared to the ones of 4MD is 

J^h( v ) d tvdx < e \\dtv\l + \ll{v)\ 2 2 . 

But, thanks to ©, we know that 

\\ X h( v )\\l < 2\v-ll{v )\ 2 + 2\u\l <2 c 0 h 2 \\Vu \\ 2 2 + 2ci h 4 \\Au \\ 2 2 + || u\\. 


Therefore, we deduce that ||I^(i >)||2 is integrable in time. 


6.2. Continuous dependence on the initial data and uniqueness of solu¬ 
tions. The last point to make clear is the uniqueness and the continuous depen¬ 
dence on the initial data for system 1 S 2 . This can be shown using either S w or 1 S 2 . 
We use the latter system. Let u be a strong solution of system and V\ and V 2 

two solutions of £2 and w = v 1 — i >2 • In order to prove the continuous dependence 
of the solutions on the initial data, and therefore the uniqueness, we proceed as in 
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section F5T21 Clearly, w satisfies (EH) and (l62l) with I 2 instead of I\, and Vi and a 2 
instead of u and v respectively in (1621) . Therefore, we have clearly 


1 d_ 

2 dt 
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~T~ 
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1 -1 
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[ (a Ko y \ 

- H 


Ms 


II Via 1 2 
(87) 


Now, on the one hand, for all 1 < a < |, using (l63l [Ml (with e = 1), [671 EE) and 
EH, we have 
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( 88 ) 


On the other hand, for all | < a < 2 , using (1631 [Ml (with e = 1 ), [671 EE (with 
e = 1 )) and (1751) . we obtain 
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(89) 


Therefore, adding (1571) to (1551) and (1571) to EH- and optimizing in the e's, we obtain 
using the fact that 15/xmax{co, ^/cj"} h 2 < v 


d_ 

dt 


ia(i)|ffi 


+ v || Aia ||2 + |(|i’i| + \v2\) a w\\l 

< [5 3 (t)-/i + max{5 2 (t),5i(t)}] ||ia(t)||^i, 









































32 


P. A. MARKOWICH, E.S. TITI, AND S. TRABELSI 


where for all 1 < a < |, we set 
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By virtue of Gronwall’s inequality, we can write 

\\w(t)f H i < \w(t = 0 )||^ e /o[^( S )-M+max{5 2 ( s ),5 3 ( S )}]d S _ 


(90) 


Next, recall that ||ui(f)|| ff i and ||'C 2 (t)||#i are uniformly bounded. Furthermore, 
integrating (1521) or (1551) with respect to time, it is rather easy to see that 

r T 

/ |Au>(s)|| 2 ds < +oo. 


In particular, we infer that 

r T 


[ ||Aui(s )|| 2 ds < +oo and f |Au 2 (s )|| 2 ds < +oo. 

Jo Jo 

The continuous dependence on the initial data and the uniqueness follow from (1901) . 
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